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We fit the ultra high energy cosmic ray spectra above 10 19 eV with different injection spectra at 
cosmic ray sources that are uniformly and homogeneously distributed in the Universe. We conclude 
that the current UHE spectra are consistent with power laws of index a between 2.4 and 2.7. There 
is a slow dependence of these indices on the cosmological evolution of the cosmic ray sources, which 
in this model determines the end of the galactic cosmic rays spectrum. 
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I. INTRODUCTION 



The current results on the energy spectrum of the high- 
est energy cosmic rays are not fully consistent. The two 
high statistics experiments, AGASA Q and HiRes Q, 
do not agree on the normalization of the ultra high en- 
ergy cosmic ray (UHECR) spectrum. In addition, HiRcs 
results are consistent with a GZK |3( suppression, and 
AGASA claims a spectrum extended to higher energy. 
With the current statistics the differences are not very 
significant - the number of events above 10 20 eV differs 
by less than 3a Q . The normalizations of the spectra are 
quite different, especially in the common E 3 x dN/dE 
presentation, but a renormalization of the energy assign- 
ment by 15-20%, which is within the reported systematic 
uncertainty of the energy assignment, of both data sets 
leads to a good agreement 0, H| . The high energy exten- 
sion and exact normalization of the UHECR spectrum 
are thus not well known, but after the renormalization 
both experiments show the same spectral shape between 
10 18 - 5 and 10 20 eV. 



There have been recently several attempts 
Hol fTTl Il2l Il3| to explain this shape with different injec- 
tion spectra of extragalactic protons after propagation to 
the observer from isotropically and homogeneously dis- 
tributed sources. The assumed injection spectra and to 
certain extent the cosmological evolution of the sources 
determine the shape of the extragalactic cosmic ray spec- 
trum at Earth. A subtraction from the observed cosmic 
ray spectrum in these models determines also the end of 
the galactic cosmic ray spectrum. There are two types 
of solutions. Flat injection spectra, dN/dE — AE~ 2 , 
are suggested in Refs. 0, 0- In this case the galactic 
cosmic rays spectrum extends above 10 19 eV. The other 
popular solution is to use much steeper injection spec- 
tra with spectral indices a = 2.6-2.7. Such solutions set 
the end of the galactic cosmic ray spectrum at lower en- 
ergy. The assumption that the extragalactic cosmic rays 
have a large fraction of nuclei heavier than Hydrogen |l4| 
modifies the propagation process and can potentially ini- 
tiate a new class of models, in which the composition at 
injection introduces more model parameters. 

Cosmic ray data is not yet good enough to prove that 
any of these models are correct. Fitting the spectra with 



several different parameters and the uncertain knowledge 
of the cosmic ray composition at 10 19 eV makes all mod- 
els plausible. 

We describe an attempt to fit the renormalized 
AGASA and HiRes spectra with injection spectra and 
cosmological evolution parameters covering practically 
the whole phase space and to use the quality of the fits 
as a measure of the plausibility of the models. This in- 
volves several assumptions that, although used in many 
previous publications, may not be correct. The main 
one is that extragalactic cosmic rays are protons and 
that the renormalized experimental spectra represent 
correctly the shape of the UHECR spectrum. One of 
the fit parameters is the required extragalactic cosmic 
ray luminosity at present time. The uncertainty of the 
luminosity depends on the arbitrary renormalization of 
the experimental data, as well on the cross correlation 
with other parameters and assumptions. 

This paper is organized as follows: Section 2 describes 
the fitting procedures that we use, Section 3 gives the 
main results of the fits, and Section 4 discusses the re- 
sults, compares to other fits and derives the main con- 
clusions from this research. 



II. FITTING THE UHECR SPECTRUM 

It was shown in Refs. 0, 0| that a renormalization 
of about 15% of the energy assignment of the AGASA 
and HiRes events, would bring the two spectra in very 
good agreement in the energy region below 10 20 eV in 
both normalization and shape. In Ref. Q it was also 
shown that the statistics of events above 10 20 eV is too 
small to achieve a conclusive result about the end of 
the UHECR spectrum. In this paper we study how 
the best fit to the injection spectrum depends on the 
source parameters: injection spectrum, luminosity and 
luminosity evolution with redshift. To do so, supported 
by the above-mentioned findings, we fit these injection 
parameters to the energy-shifted spectra in the energy 
range 10 19 -10 20 eV. We chose as lower bound 10 19 eV be- 
cause we expect at this energy the particles to be mostly 
extra-galactic and 10 20 eV as higher bound because of 
the sparseness of data above this threshold. 

To calculate the expected spectrum from a isotropic 
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homogeneous distribution of proton sources we use the 
analytical approach presented in Ref. 0, |r|- In this 
approach all proton energy losses are included (redshift 
losses, pair and pion production losses) and they are all 
treated as continuous. This is the correct treatment for 
redshift losses and it is well suited for pair production and 
for pion production at large propagation distances. At 
small propagation distance, however, the large inelastic- 
ity of the pion production process produces large fluctu- 
ations in the expected fluxes that cannot be reproduced 
in the continuous energy loss approximation. A Monte 
Carlo simulation is better suited in this case. In the 
present paper, since pion production affects mostly the 
highest part of the energy spectrum, and we are only in- 
terested in the spectra below 10 20 eV, we can safely use 
the continuous energy loss approximation for all energy 
loss processes. 
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FIG. 1: Solid line: loss length for photo-pion and photo-pair 
production for protons as calculated in Ref. 0. The dashed 
lines report the separate contribution of the two processes. 
The dotted line shows the loss length for redshift losses. The 
squares represent the loss length of Ref. |2Cj|. 

We assume the sources inject protons with a power-law 
spectrum, E~ a , with a sharp cutoff at E max = 10 21 - 5 eV. 
Changing the value of -Emax does not appreciably af- 
fect the results in the energy region we are interested 
in. We assume the source emissivity to evolve as L(z) — 
Lq(1 + z) m , where Lq is the present cosmic ray emissiv- 
ity of the sources and m = corresponds to absence of 
evolution. There is no cutoff to the evolution of the lumi- 
nosity because we are only interested in the region above 
10 18 5 eV and in this region the contributions to the ob- 
served flux can come only from sources with z < 0.5. 
To calculate the energy losses we use the loss-lengths of 
Ref. U, plotted in Fig. □ that were shown to be indis- 
tinguishable from the ones of Ref. [2(j in the low en- 
ergy region and within 15% at high energy |(|. We as- 
sume a ACDM universe with £7a = 0.7, £l m = 0.3 and 



H a = TSkms^Mpc- 1 . 

We explore the parameter region in a from 2.05 to 3.00 
in steps of 0.05 and in m from to 4 in steps of 0.25. For 
each (a, m) pair we calculate the expected flux and then 
the best fit emissivity, Lq, minimizing the \ 2 indicator. 
To emulate the experimental energy resolution we include 
30% Gaussian error distribution in the spectra after prop- 
agation. The effect of the energy resolution is to smooth 
the features produced by the propagation on the photon 
background and to lessen the GZK suppression [4|. 
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FIG. 2: Upper panel: Best fits to the AGASA -15% dataset 
in the 10 19 -10 20 eV energy range. Squares with error bars: 
HiRes dataset. Circles with error bars: AGASA&Akeno 
dataset. Filled symbols: points used in the fit. Smaller data 
points: data from HiRes2 and from Akeno (not shifted) to 
show the low-energy shape of the spectrum. Solid lines: best 
fits for different values of m. Dotted lines: galactic compo- 
nent needed in order to fit the spectrum at lower energy. The 
numbers attached to the lines indicate the value of m. Lower 
panel: best fit for the HiRes +15% dataset. 

For AGASA and Akeno we use the data from Ref. pj, 
whereas for HiRes 1&2 we use the data from their web- 
site |U which is very close to the published results 
Following the suggestion of Ref. we shift the AGASA 
and HiRes energies respectively by —15% and +15%, 
while we leave the Akeno energies unchanged j^. To 
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do the shift we proceed in the following way: we calcu- 
late d i g^ (£ ) (E) and we assign this flux, calculated in 

E, to the energy kE. This means that: (^(fci?)) ncw = 

dN (E))° W or (k 3 E 3 ^(kE)) ncw = k 2 (E 3 ^(E))° U , 
1.15 for HiRes 1&2 and k = 0.85 for AGASA. 
These new, shifted, datasets agree quite well almost over 
the whole energy range as shown in Fig. [21 



III. RESULTS FROM THE FITS 
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We applied the method presented in the previous para- 
graph to the AGASA —15% data fitting the points in 
the energy range 10 19 -10 20 eV. The results are shown 
in the upper panel of Fig. |2 where we plot the best fits 
for m = 0, 3, 4 (solid lines). We only show results for m 
= 3 and 4 because these values bracket the cosmologi- 
cal evolution derived from star forming regions and from 
gamma ray bursts. The corresponding slopes are, respec- 
tively, a = 2.55, 2.45, 2.45 with la errors of about 0.20. 
As it is clear from the plot all the three curves fit well 
the data in the region considered, with different degrees 
of goodness at low energy. The dotted lines represent 
the needed galactic component to fit the spectrum. The 
best fit with m=4 does not allow for galactic component 
above 10 17 eV. 

We repeated the same exercise for the HiRes +15% 
dataset (which is shown in the lower panel of Fig. [2J 
and the results were similar. The slopes are somewhat 
steeper, by about 0.05-0.1, well within the similarly large 
uncertainties. The main difference is at energies much 
lower than the fitting range, where the fits of the HiRes 
data set do not allow for a galactic component in cases 
with cosmological evolution. The flux of extragalactic 
cosmic rays below 10 19 eV has to be slightly decreased 
by some additional process in order not to exceed the 
Akeno and HiRes measurements. The shape of the spec- 
tra in the considered region is, however, the same for 
both experiments. 

It has to be noted that the inclusion of the error dis- 
tribution in the fit affects the spectral shape - the pile-up 
approaching 10 20 eV is visibly smoother. Since the points 
immediately above 10 19 eV have the lowest error bars, 
and thus affect the fit the most, the slope of the spec- 
trum is increased by at most 0.05, much smaller than the 
uncertainties from the fits. 

The best fitting parameters depend slowly on the fitted 
energy range. Fig.Elshows the fit of the HiRes +15% data 
set for the energy range between 10 185 -10 20 e V range. 
The best fit with m=4 again does not leave space for 
a galactic component below 10 18 eV. The best spectral 
indices are 2.6, 2.5, and 2.5 respectively for m = 0, 3, 
and 4, slightly smaller than those for the higher fitting 
threshold and almost identical to the AGASA —15% set 
shown in Fig. [3 The la error bars decrease to about 
±0.1. The effect of the wider fitting range on the AGASA 
— 15% spectrum is similarly small, although for this data 
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FIG. 3: Best fits to the HiRes +15% dataset in the 10 18 ' 5 - 
10 20 eV energy range. Symbols are the same as in Fig. 
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FIG. 4: Contour plots of log 10 \ 2 as a function of a and m for 
the AGASA —15% dataset (upper panel) and for the HiRes 
+15% data set (lower panel) fits above 10 19 eV. The white 
lines are the contours corresponding to la, 2a and 3a. 

set it increases the spectral slopes for all m values by 
0.05-0.1. The la fit errors decrease to slightly less than 
±0.1. 

As a consistency check we also fitted the unmodified 
AGASA and HiRes spectra. The results we obtain are 
much like the ones presented above. The best fit pa- 
rameters differ by about 0.05-0.1, with the same ±0.2 
errorbars. 

We performed several other fits varying the fitting 
threshold between 10 185 and 10 19 2 eV and convinced 
ourselves that all fits returned consistent results within 
the la errors of the presented fits as shown in Fig. 0] 

In this figure we plot in the top panel log 10 x 2 &s a func- 
tion of (a,m) for the AGASA -15% fit above 10 19 eV, 
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FIG. 5: Required emissivity as a function of (a,m). Lines: 
best fit emissivity (above 10 19 eV) as a function of a. The 
numbers attached to the lines give the value of m. Every line 
is highlighted with different thicknesses corresponding to the 
confidence bands for la, 2a and 3a. The black squares are 
the best fits obtained in other works. W&B corresponds to 
Ref. HEI, BGG corresponds to Ref. @. 

i.e. for six degrees of freedom. The analogous fitting of 
the HiRes +15% set is shown in the bottom panel. The 
white contours are the confidence bands for la, 2a and 
3a. All to values can provide a good fit to the data as the 
best fit a value slowly decreases with increasing to. This 
correlation is easily understood as with a smaller value of 
a less low energy particles are injected and to compen- 
sate for that one needs a stronger evolution of the sources 
to increase the number of low energy particles reaching 
the observer. For the HiRes +15% dataset the best fit 
parameters are in the strip connecting (a = 2.6, m = 0) 
and (a = 2.4, to = 4). As it is clear from the confidence 
bands in the plot, the present data sets do not restrict 
very much the values of the parameters, a being deter- 
mined with an uncertainty of ±0.2 for a given to and to 
being almost free for a given a. It is obvious, though, 
that fits with a flat injection spectrum do not give good 
X 2 values even with a strong cosmological evolution. In- 
jection spectrum with a = 2.0 would be in the 3a range 
only in the case of m — 4. Flat injection spectrum mod- 
els require that the galactic cosmic ray spectrum extends 
tol0 195 eV. 

If the shape of the cosmic ray spectrum is the same as 
the one derived from the existing experimental statistics, 
even much higher future statistics from the Auger obser- 
vatory 23] would not help to solve it. We performed a 
fit with the current spectral shape and increased statis- 
tics that corresponds to the one expected from Auger. 
The la errors on a for the fits above 10 19 eV became 
only slightly narrower ±0.15. The measurement of the 
cosmic ray chemical composition, or, a measurement of 
the flux of cosmogenic neutrinos generated by UHECR 
in propagation to us 01 > are needed to disentangle the 
two parameters. 

In Fig. [S] we plot the best fit present day emissivities 
above 10 19 eV as a function of a for different values of 
to. In this plot we also show the values obtained in 
Ref. USUI. The differences with the results of Ref. 6] 
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FIG. 6: One a error band (shaded) of the AGASA fit with 
m=0 above 10 19 eV. 

are likely due to the slightly different dataset, to the dif- 
ferent range of data used for the fits and to the inclusion 
in our calculation of the experimental energy resolution. 
There is also a factor because the fits were performed 
with a —10% shift of the AGASA data set instead of the 
— 15% shift used here. It is interesting to note that the 
required values of the emissivity above 10 19 eV cover a 
narrow range between 2 and 6 • 10 44 erg Mpc -3 yr _1 and 
that the required luminosity increases with the flattening 
of the injection spectrum. This is a consequence of the 
fact that we only present the luminosity required above 
10 19 eV. If we were to extend the energy spectrum to 
lower energy, say to 10 17 eV, we would observe exactly 
the opposite trend - steeper injection spectra would re- 
quire much higher luminosity than flatter ones. 



IV. DISCUSSION AND CONCLUSIONS 

After fitting the shifted AGASA and HiRes data sets 
in terms of injection spectral index and cosmological evo- 
lution of the cosmic ray sources for an isotropic and ho- 
mogeneous source distribution we obtained current emis- 
sivities above 10 19 eV that differ only by about a factor 
of two. In this sense we confirm the statement of Wax- 
man flfil ] that approximately the same emissivity is re- 
quired for a wide range on injection spectral indices. We 
disagree with the estimate of the central spectral index 
m Ref. 

and find a significantly steeper one. 
Best fit spectral indices are, however, not well re- 
stricted by current statistics. In Fig. H3 we show with 
shaded area the la errors of the best fit prediction from 
the AGASA data above 10 19 eV (top panel of Fig. 0) 
for 77i = 0. The figure emphasizes the perils of all fits 
of the extragalactic cosmic ray component with the cur- 
rent statistics. Such fits are most sensitive to, and at- 
tracted by, a small number of experimental points with 
the best statistics, in our case four points between 10 19 
and 10 19 4 eV. Fitting uncertainties do not affect much 
the higher energy spectra where the GZK suppression 
prevails almost independently of the injection spectral in- 
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dex but create a large uncertainty below 10 185 eV. This 
uncertainty makes estimates of the end of the galactic 
cosmic ray spectrum by subtraction of the model predic- 
tions from the total observed flux unreliable. 

The luminosities that we show in Fig. [5] apply only to 
energies above 10 19 eV. If one is interested in the to- 
tal cosmic ray luminosity of the cosmic ray sources one 
should continue the integration to much lower energies. 
This introduces several possible new astrophysical pa- 
rameters that come from the exact acceleration mecha- 
nism of the extragalactic cosmic rays. One could inte- 
grate down to the proton mass and obtain the highest 
possible emissivity. On the other hand studies of parti- 
cle acceleration at relativistic shocks |24j find a minimum 
acceleration energy of mFg hock which for r s h oc k = 1000, 
as in gamma ray bursts, could be 10 15 eV and would de- 
crease significantly the required emissivity. Modifications 
of the cosmic ra y sp ectrum on propagation, such as sug- 
gested in Refs. [l(J, fill l2q| because of magnetic horizon 
of the high energy cosmic rays [2(| would not change the 
required emissivity. Such modifications would however 
change very much the extra-galactic cosmic ray spectrum 
suppressing the flux at the lower energy end and by con- 
sequence changing the shape of the end of the galactic 
cosmic ray spectrum required in order to fit the observa- 
tions. 

Fits of the monocular HiRcs data have been performed 
by the HiRes group . The best fit is obtained for a — 
2.38±0.04 and m = 2.8±0.3 for a total of 42 data points 
above 10 17 eV and respectively 39 degrees of freedom. In 
additon to the different energy range of the fit, HiRes 
assumes a 'toy' galactic cosmic ray model based on their 



composition measurement |l9j that suggests domination 
of the extragalactic cosmic rays above 10 18 eV. 

The HiRes fit is probably dominated by lower energy 
cosmic rays 

( 10 17.5-18.5 e y) with much sma n er 

error bars. 

The two fits give different central values for a and m 
but they are qualitatively consistent in the conclusion 
that even with a strong cosmological evolution of the 
cosmic ray sources the observed spectra do not support 
fiat injection spectra. 

We have fitted the shape of the ultrahigh energy cosmic 
ray spectrum above 10 19 eV assuming that these cosmic 
rays are protons, and that the sources of these protons 
are uniformly and homogeneously distributed in the Uni- 
verse. The fits of the scaled AGASA and HiRes data 
sets allow for power law injection spectra in the range 
AE~( 2A ~ 2 - 7 ^ for cosmological evolution of the cosmic ray 
sources between (1 + zf to (1 + z)°. The cosmic ray 
emissivities above 10 19 eV required by different models 
are within about a factor of two in this range. The best 
fit spectral index decreases for strong evolution models. 
Flatter injection spectra do not fit well the cosmic ray 
spectra above 10 19 eV. This also means that the end of 
the galactic cosmic ray spectrum is at, or below, 10 18 eV 
depending on the cosmological evolution of the extra- 
galactic cosmic ray sources. Consistent data on the cos- 
mic ray composition in the energy range above 10 17 eV 
are required in order to reveal the end of the galactic 
cosmic ray spectrum and thus help determine that of the 
extragalactic sources. 
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